Statics and Dynamics of Condensed DNA within 
Phages and Globules 

Theo OdijJ* 
Theory of Complex Fluids 
Faculty of Applied Sciences 
Delft University of Technology 
Delft, the Netherlands 

Abstract 

Several controversial issues concerning the packing of linear DNA in bacteriophages and glob- 
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uniform density. A new electrostatic model is introduced for computing the osmotic pressure of 
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Coulomb sliding friction is proposed. A general discussion is given of our current understanding of 
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cyclization within globules and single-molecule determination of the loading forces. 



*E-mail: odijktcf@wanadoo.nl 
Address for correspondence: 

T. Odijk, P.O. Box 11036, 2301 EA Leiden, the Netherlands 

1 



1. Introduction 

The compaction of DNA in biological cells and viruses presents us with rather unconven- 
tional problems in the physics of soft matter. In the case of phages or viruses, the DNA may 
be so close-packed that the thermodynamics is no longer extensive because it is dominated 
by the energy arising from regions of tight bending (Odijk 1998). Thus we are dealing with 
a defect rather than a bulk system. Confronting the mathematical physics of the DNA con- 
figurations as such is a formidable undertaking. Still, the condensed DNA globule ought to 
have relatively minor fluctuations in the density (Ubbink & Odijk 1996) so that simplified 
analyses are possible. 

A general maxim in theoretical biology is not to theorize before one has assimilated per- 
tinent data and the same would seem to apply to the field of DNA compaction also. In 
setting up the computations outlined here, I have relied heavily on what we have learned 
from a number of recent experiments. The spoollike structure of DNA in T7 and T4 bac- 
teriophages has been established in a quantitative fashion (Cerritelli et al 1997; Olson et 
al 2001). For the phi29 phage we now have a fairly good comprehension of the connector 
motor (Simpson et al 2000) and the typical forces involved in DNA packaging as a function 
of length reeled in (Smith et al 2001). Cyclization experiments (Jary & Sikorav 1999) of 
DNA in condensed globules allow us to estimate the local friction the DNA coil experiences 
in a congested environment. 

Provisional theories for DNA compacted in phages were already proposed a couple of 
decades ago (Reimer & Bloomfield 1978; Gabashvili et al 1991; Gabashvili & Grosberg 
1992). Nevertheless, the notion that the tight curvature of the DNA sheath surrounding a 
hole determines the interaxial spacing, was stressed only recently (Odijk 1998). Theoretical 
activity in this area is currently gaining momentum but will be discussed in the last section. 
The experimental literature on DNA packaging in phages and viruses is vast and compli- 
cated and, owing to space limitations, cannot be dealt with here (for recent reviews of the 
experimental situation, see Grimes et al (2002) and Jardine & Anderson (2003)). Here I 
wish to emphasize several theoretical issues that have not been focussed on before. 

A continuum approximation of the curvature energy of DNA condensed within a capsid 
is derived in terms of the inhomogeneous DNA density in the limit of small fluctuations. 
A further useful approximation for close packed DNA is to let the density be constant. 
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It turns out that it is then possible to derive general nanothermodynamic relations for 
quantities like the osmotic pressure, force of insertion, etc. Owing to the very high degree 
of packing in phages, we are obliged to reinvestigate the classical cell model for rodlike 
polyelectrolytes. The interpretation of the cyclization experiments of linear globular DNA 
leads to the (surprising?) conclusion that simple reptation seems to be valid though I note 
that regions of tight curvature are often absent in globules. The friction of DNA in highly 
congested phages is discussed in terms of current thinking on sliding friction at nanoscales. 
In the last section, I confront these new computations as well as previous theoretical work 
with the experiments mentioned above. 

2. Continuum approximations 

The persistence length of DNA is generally of the order of magnitude of the typical size of 
bacteriophages. This makes the statistical treatment of DNA viewed as a confined wormlike 
chain with self- interaction very difficult. Fortunately, there is one limit where we can make 
considerable headway, viz. when the DNA is close-packed at moderate to high densities 
owing to bending stresses. Such a condensed state may exhibit so-called frustration because 
the packing - hexagonal for instance - cannot be perfect (Earnshaw & Harrison 1977; Pereira 
& Williams 2000), but this effect will be neglected here. 

Let us then consider a configuration of a DNA molecule enclosed within a capsid of 
volume Vq and defined by the radius vector f(s) as a function of the contour distance s from 
one end. The bending energy U c in a Hookean approximation is written in terms of the local 
radius of curvature R c (s) of the DNA curve at s (Yamakawa 1971) 



Here, the bending force constant equals PksT where P is the DNA persistence length, k B 
is Boltzmann's constant and T is the temperature. We wish to introduce the normalized 
density of DNA segments at position f expressed in terms of a delta function 




(1) 




(2) 
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The contour length L of the DNA consists of N "segments" of length A, the distance between 
elementary charges viewed along the DNA helical axis. We next insert the density into eq 
(1) and average over all configurations (denoted by ( )) which is possible because f(s) is a 
unique nonintersecting curve. 



We have also introduced a close-packing approximation (p c (r) / 'R 2 {r)) ~ (p c )(R~ 2 ) = 
p{r)R~ 2 {r) which is certainly accurate when the DNA undulations away from the average 
are small. I note that eq (3) pertains to both purely mechanical and statistical mechanical 
types of theories. Even in the former, we still need to average over all configurations, for the 
capsid has a complicated shape so there must be a host of DNA conformations of virtually 
identical energies. The bending energy given by eq (3) was used recently (Odijk and Slok 
2003) to set up a simple density functional theory of packing and will be commented on 
below. 

Next, we consider a more severe approximation by letting the density be constant. It 
is conveniently expressed in terms of the area S of the unit cell of the packaged DNA 
(p = l/SA) which we often take to be hexagonal. The supposition of uniform density 
allows us to deduce expressions of some generality. The bending energy is rewritten as 



where the infinitesimal volume dr = Sds. The volumes Vj and Vq are the inner and outer 
volumes, respectively, enclosing the packed DNA (see fig 1). We assume there is only one 
hole though it is stressed that general shapes of the capsid and hole are allowed. Now it can 
be argued by scaling and other arguments (Odijk 1986; Selinger & Bruinsma 1991; Kamien 
et al 1992) that the twist and splay contributions to the free energy of a tightly packed 
liquid crystal are small compared with the bending energy. Thus we interpret R~ 2 as an 
energy density of pure bending which depends on the director n{f) but not on S, the size 
of the unit cell. Moreover, because there should be no hairpins at close packing, the DNA 
mesophase is splayless at constant density (de Gennes 1977). 




(3) 





(4) 
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There is also a free energy of the DNA interacting with itself. We omit minor curvature 
and surface contributions (Ubbink & Odijk 1995) and express this energy as an extensive 
form proportional to the contour length L 



F mt = Lf(S) 



(5) 



This pertains to a hypothetically straightened DNA lattice. Accordingly, we now want to 
minimize the total free energy 

F tot = U + F int (6) 

subject to a volume constraint 

Vi = V - LS (7) 

We minimize F tot with respect to the director n (implicit in Rr 2 ) and to the cell size S, 
keeping the outer volume Vq fixed but also the chemical potentials of the small ions of the 
buffer in the reservoir containing the phage. 

SF tot 



5n 
dF tot 







(8) 



S,Vi(S) 



OS 



These expressions yield 



dl(fi 



dS 







(9) 



in terms of an optimized inverse square radius of curvature _R m 2 and an integral over the 

— * 

surface of the inner hole (see fig 2). A variation in S induces a change in the vector I at 
point fj, pointing normal to the inner surface in the direction of the hole. The variation of 
the inner volume can be similarly written with the help of eq (7) 



dS 



-I 

J A, 



dfi 



dl{n 



OS 



-L 



(10) 



Hence, it is expedient to define a measure of the hole size given by 



E7 2 = 





n) 


dl(n) 
dS 


4 df * 


dl(fi) 
dS 





(11) 
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Noting that —f'(S) is simply the osmotic pressure n os , we finally attain the nanothermody- 
namic relation 

n ° s = -SL + ^SEf (12) 
The usual thermodynamics does not hold: if the hole is small, we have Vo = O(SL) and 
U = 0(PV^ 3 k B T/S) and so the osmotic pressure is not a purely intensive quantity. 

Next, the average pressure V on the capsid wall exerted by the DNA is computed by 
considering a virtual change in Vq. Concomitantly, the director configuration n(r), the inner 
volume Vi(Vo, S) and cell size S also change in such a way that the free energy remains a 
minimum (eq 8) and the constraint (eq 7) is obeyed. The resulting analysis is similar to 
that given above, leading to 



v = dF tot dU 



Pk B T Pk B T U Pk B T 



s 2SE? 2SE 2 SL os 2SE 2 



+ *os ~ ITW^o ( 13 ) 



dV dV 

The typical size E of the capsid is defined via its surface similarly to eq (11). I note that this 
pressure is not at all identical to the osmotic pressure as one might have naively surmised. 

It is also interesting to study a quasistatic force T (in the absence of dissipative losses; 
Kindt et al 2001) upon increasing the length of DNA in the capsid. We require local 
equilibrium (eqs (14)) as n, S and Vi(S, L) adjust to the changing degree of packing subject 
to eq (13). It is straightforward to show that the "pressure" T jS of insertion is given by 

l-^l dF tot f U 

-T = ML = S +7Tos+ SL (M) 

Again, it would have been difficult to anticipate this expression. We have derived eqs (12- 
14) for a phage of general shape at the expense of introducing the capsid and hole sizes as 
unconventional averages. 



3. Hexagonal DNA lattice at high concentrations 

We have argued that we merely need the free energy per unit length f(S) of a straightened 
DNA lattice in order to evaluate pertinent quantities of the DNA enclosed in a capsid. In 
this section we assume such a lattice is perfectly hexagonal without defects. Actually, 
the DNA chains also undulate about a reference configuration but this effect is temporarily 
disregarded. The lattice model for rodlike polyelectrolytes is, of course, well known (Oosawa 
1971), but the usual cell model which disposes of the hexagonal symmetry, breaks down 
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at high degrees of packing when the polyion cylinders are close. A new approximation 
is introduced here, though some general features of the Donnan equilibrium need to be 
discussed first. 

For now, the lattice is free from simple electrolyte and consists of parallel cylindrical 
polyions together with attendant counterions dissolved in the water in the intervening void. 
In some model or approximation one introduces, there is generally an interface or line of 
symmetry where the negative electrostatic potential is a maximum and is set equal to zero. 
The electric field vanishes so the osmotic pressure n os and chemical potential of the ions 
reduce to (Israelachvili 1985) 

T^os = Pik B T 

Pi = Pref + k B T In ^ (15) 

where pi is the concentration of counterions at the interface. We restrict ourselves to monova- 
lent ions bearing an elementary charge q which are ideal save for their Coulombic interaction 
mediated by a supposedly uniform permittivity e. 

Next, the lattice is brought into thermodynamic contact with a large reservoir containing 
monovalent salt of concentration c s . The activity coefficients of the small ions are set equal 
to unity so the resulting Donnan equilibrium, due to the equality of the chemical potentials 
of the small ions in the respective phases, yields 

p_(p l + P+ ) = c 2 s (16) 

The polyelectrolyte suspension contains positive and negative ions arising from the simple 
salt, whose densities must be set equal in view of electro neutrality (p+ = p). Therefore, 
upon eliminating p + and p_, the osmotic pressure is simply (Oosawa 1971) 

n s = {pi + P+ + P- - 2c s )k B T 

= p l k B Tg{w) (17) 

9 (y) EE (l + y*)V*-y 

w = 2c s /pi 

In the so-called cell model (developed by Fuoss, Katchalsky, Oosawa and others, see Oosawa 
1971), the counterion concentration pi at the interface in the salt-free case is given by 
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This is valid at low and intermediate concentrations. For DNA, we have a segment length 
A = 0.17 nm and the Bjerrum length Q = q 2 / eksT = 0.71 nm at room temperature implying 
that a fraction of counterions appears to be associated with the DNA rods. This is often 
termed counterion condensation (Manning 1969). 

We now address the Poisson-Boltzmann equation in the lattice at high concentrations (in 
the salt-free case). The ion distribution p(r) is connected to the (negative) electric potential 
ip via a Boltzmann distribution 

Pi = qpiexp(-q(p/k B T) (19) 

which is inserted in the Poisson equation 

Atp = (20) 

This procedure is not entirely rigorous although it is quite accurate for monovalent ions (the 
difficulty is that the two densities in eqs (19) and (20) are very similar though not identical; 
see Fixman 1979). 

Let us consider the case when the surfaces of the polyion cylinders are separated by only 
a short distance h which is much smaller than the cylinder radius a (see fig 3). Hansen et 
al (2001) then still insist on using a cell model in which the effect of the rods surrounding 
a test rod is mimicked by a cylindrical boundary. This seems a rather severe approxima- 
tion. An alternative is to view the hole enclosed within the triangle in fig 3 as a separate 
entity which is feasible provided the electrostatic screening is sufficiently high. Solving the 
Poisson-Boltzmann equation in this geometry is still difficult, so we replace the hole by three 
sections involving thin layers of electrolyte, and one effective hollow cylinder whose radius b 
is appropriately chosen. I consider the latter first. The boundary condition for the electric 
field at its surface is 

dip 



dr 



(21) 



e 

with r the radial coordinate. The negative charge density a q = —qa with a the number 
of charges per unit area to be chosen below. We have to solve eq (19) and (20) so it is 
convenient to introduce the dimensionless variables \I/ = qp/ksT and f = r/A where the 
screening length A is analogous to the Debye length 



A 2 = 4nQpi 
8 



(22) 



The counterion density pi now pertains to the centreline of the cylinder. There, the potential 
is set equal to zero and the electric field is also zero so that eqs (15) hold. It is then quadrature 
to solve the Poisson-Boltzmann equation: 

^ = 21n(l-if 2 ) (23) 



The concentration is determined from eqs (21) and (22). 



b\ 2 8ttA 



A / 1 + ttA 
A = Qab 

Hence, the osmotic pressure is conveniently expressed as 

7T os = pik B T 

1ok B T 



(24) 



(25) 



1 + TTA, 

At small A, this reduces to a condition of electroneutrality as it should, since no counterions 
are condensed in this limit. 

To progress, we need to set the variables a and b. The size of the thin boundary layers 
discussed below is so small that they may be neglected in the following argumentation. The 
triangular region in fig 3 bears half an electron charge per length A along the major axis. 
The surface density of the effective cylinder is thus a = 1/AnbA implying that A ~ 2 is 
independent of the interaxial spacing H (Q/A = 4.2 for DNA). The radius b is chosen such 
that the volumes of enclosed liquid are the same in the respective compartments. 

irb 2 = \V3H 2 - -ira 2 (26) 
4 2 

We rewrite eq (25) by introducing the volume fraction v = Po/pdna of DNA (the segment 
density of pure DNA = Pdna] the area of a hexagonal unit cell is S — \/3H 2 /2). 



k os = PiksT 
Po k B T 



2(1 -v) 

The density p is identical to the uniform DNA density discussed in the previous section 



(27) 
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Obviously, our choice of cylinder parameters is not at all rigorous so the numerical coef- 
ficient in eq (27) is tentative. The scaling of parameters in eq (25) is convincing, however, 
because the correct limit given by eq (18) is attained at low densities (H 3> a i.e. A » 1). 
Furthermore, the screening length A is indeed small. For instance, at a separation h very 
close to zero, we have b = 0.23a from eq (26) and since pi ~ l/47ib 2 A we get A ~ 0.47a from 
eq (24). Hence, the respective triangular sections (fig 3) in the entire lattice basically act as 
independent units as far as the electrostatics is concerned. 

The thin boundary layers are treated in a Derjaguin approximation taking advantage of 
the relatively large radius (a > A). Of course, discrete charge effects play a role now, so our 
analysis is not quantitatively accurate. For two flat plates of charge density o = l/2naA 
equal to that of DNA, we have 

Kos = — ^ — (28) 

The scaled density A p = aQh, the analogue of A in eq (24), is much smaller than unity 
so the counterions are essentially uncondensed. The strategy of Derjaguin is to superpose 
normal forces exerted by opposing infinitesimal platelets on the two nearby charged objects 
(Israelachvili 1985). We therefore focus on two particular cylinders in fig 3 and introduce 
the plane defined by their centre axes. Two interacting platelets are located at a height x 
from this plane and the distance Z between them is Z = h + 2z with z = a — (a 2 — x 2 ) 1 / 2 . 
Hence, for relatively large a, we have Z ~ h + (x 2 /a). The normal force between the two 
platelets per unit area is from eq (28) 

T r (Z) = 2 -^f (29) 

and the total force between the two cylinders of length / is hence given by 

/•oo 

T(K) ~ 21 / dx T V {Z) 

J Z=h 

Aalk B T f 00 , / a; 2 V 1 
da; 1 + 



h In \ ha 



= 2naa 1/2 h- 1/2 lk B T (30) 

This can be integrated to give F 2 , the free energy of interaction between the two cylinders. 
If there are n rods in the actual hexagonal lattice, the total free energy is F int = 3nF 2 for 
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there are six cylinders surrounding a test rod though we have to avoid doublecounting. By 
definition the osmotic pressure is 

1 dF mt 



3 1 / 2 nLH dH 



AHa 1 / 2 !! 1 / 2 
3Hp k B T 



(31) 



2a 1 /2/ l i/2 

upon using the density for DNA, a = Xjl-naA and eq (30). 

At small separations the density pi in the middle of the hollow cylinder (eq (27)) is 
the maximum density within the triangle in fig. 3. There is another p i 2 defined midway 
between the surfaces of two nearby cylinders. In accordance with eq (28), this is given by 
p~i 2 ~ 1/nahA so the screening length A = Aha/AQ from eq (22) is quite small, significantly 
smaller, in fact, than the transverse scale a l l 2 h l l 2 of the "surface of interaction" evident in 
the integrand of eq (30). This justifies the Derjaguin approximation but also the splitting 
up of the original triangular region into a hollow cylinder and three thin boundary layers. 
We now place the lattice in a reservoir containing monovalent electrolyte of concentration 
c s . A Donnan equilibrium is established as described by eqs (16) and (17) with pi given by 
eq (27). 



4. Friction within condensed DNA 



As the phage connector forces the DNA genome into a capsid, the DNA coil is bound to 
twist in view of its helical backbone. Since the other end probably also becomes constrained 
at some stage, we expect a twist-to-writhe transition to occur in dealing with the DNA 
dynamics. At first, the DNA density is statistically isotropic more or less, although the 
density decreases towards the central region owing to the curvature stress (Odijk & Slok 
2003). As the concentration is enhanced, the excluded- volume effect, at some point, causes 
the DNA globule to become (probably) nematic in the outer confines of the capsid while 
the inner region remains isotropic (I disregard the chirality of the DNA interaction for 
simplicity). As more DNA is piled up into the phage, the entire globule could become 
nematic in a spoollike configuration. At the next stage, a hexagonal (or hexatic) phase 
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appears, again, in the outer region in coexistence with a nematic inner phase. Ultimately, at 
high degrees of packing, the entire structure becomes hexagonal apart from faults caused by 
frustration. This scenario is tentative of course. An a priori theory of globular organization 
is formidable as noted earlier. One problem is how the DNA loops are forced into ordered 
phases though we expect rippling transitions to occur (Odijk 1996). 

I here present two highly simplified notions of the DNA dynamics. At low degrees of 
congestion, the usual reptation dynamics (de Gennes 1979) may be thought to apply, where 
a test chain slithers through a tubelike environment. Gabashvili and Grosberg (1992) al- 
ready considered reptation during the ejection of DNA from phages. The friction the DNA 
experiences as it is transported through the connector, was also accounted for. Here, the 
reptation of DNA is discussed in the context of cyclization in order to gauge the accuracy 
of such a picture. At high degrees of congestion, the osmotic stress builds up because of the 
strong curvature stress. Simple reptation no longer seems likely then so I introduce a model 
based on sliding friction. 



Cyclization within globules 

The DNA is compacted into a globule by external forces (e.g. by adding inert polymer 
or multivalent ions, by confinement within a capsid in the initial stages). As argued by 
Ubbink and Odijk (1995, 1996), it is often possible to view the globular shape as being 
determined by balancing bending against surface forces while keeping the volume V t fixed. 
For an ideal toroid defined by the two radii R r and R t shown in fig 4, the surface area is 
A t = An 2 R r R t and the volume V t = SL = 2iT 2 R 2 R t is also given in terms of the (constant) 
area S of the unit cell. The two ends would never come into contact if the shape of the 
toroid were to remain the same, for the free energy of change in volume is prohibitive. We 
thus investigate fluctuations in the globular shape. I desregard free energy contributions 
arising from crossovers, knots and other defects (Park et al 1998; Arsuaga et al 2002; Hud 
& Downing 2001; Nelson 2002). 

The curvature stress is low if R t > 2R r (Ubbink & Odijk 1995; Ubbink & Odijk 1996). 
Under these conditions the bending energy (eq (1)) is approximated by 

u = M < 32 > 

There is a surface tension a scaled by k B T whatever the mechanism of compaction. The 
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surface free energy is then 

F s = aA t k B T (33) 

Since V t is constrained, we eliminate one of the radii so as to express the total free energy 
in units of ksT as 

F t = \PLR- 2 + 2 3 / 2 naR 1 t /2 V t 1/2 (34) 
Minimization of F t with respect to R t yields an optimum radius of the toroid 

R ^ = 2V 2 ™V (35) 
Therefore, the fluctuations in globular shape are governed by a Gaussian distribution 

G oc exp-^'Ai? 2 (36) 
<Ai? 2 > ee {{R t -R t , m f) = {F' t r 1 

This is derived by a straightforward Taylor expansion of eq (34) within the Boltzmann factor 
exp(-Ff). 

We next need to assess how the winding number n changes in response to alterations in 
shape. This is simply nS = nR 2 so a fluctuation in n is given by 5n = 2nR r>m 5R r / S. The 
volume constraint implies 2R rm R tm 5R c + i? 2 m (5i? t = 0. If the two ends are to meet, we 
require 5n > 1, leading to an inequality for the second moment (via eq (36)). We finally 
attain the condition L > 10tt 2 P which, remarkably, is independent of the surface tension. 
Accordingly, by this mechanism, cyclization can only occur if the DNA is quite long. 

A reptation time for cyclization is now readily written down by dimensional reasoning. It 
must be proportional to the friction exerted by the fluid (water presumably, with a viscosity 
770, since the globule is not very compact), i.e. proportional to the contour length L. For 
cyclization, one of the DNA ends has to diffuse a distance, typically about irR t at least, with 
respect to the other end. 

Tr ^(^* (37) 
k B l 

We do not address the dynamics of the two ends as they approach each other within a region 
of size ~ R r . 
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Sliding friction in highly compacted states 

Let us first consider a classic example of Coulomb friction (Persson 2000). A block of 
some material rests on a flat surface (fig 5). Owing to gravity, for instance, the block exerts 
a pressure J-s/Q on a section of the surface of area Q. If we now exert a force T t on the 
block, tangential to the surface, the blocks fails to move at first. It is only at a certain T t 
given by 

^-jyn- w (38) 

that the block is set in motion. The friction constant uj is well defined, at least within certain 
bounds on the velocity thus induced. 

In a highly congested DNA globule compressed by bending stresses, we may also discern 
a pressure ir os acting on a surface area 2y / 3 LH where 2\f?>H is the length of the perimeter 
of the unit cell assumed to be hexagonal. Hence, the minimum force needed to set the DNA 
in motion should be 

T t = 2\/?>ujHL'k os 

(39) 

2uLPk B T 
" HE? 

At very high degrees of packing, the first term in eq (12) may be neglected. The minimum 
force increases fast upon further packaging of DNA in a phage (increasing L, decreasing H 
and Ej). It has been assumed that the DNA slides homogeneously through its own tubelike 
environment. 



5. Discussion 



Ordered phases in the bulk and in phages 

It is important to summarize several results for liquid-crystalline DNA in the bulk that 
are relevant to DNA packaged within phages. The literature on DNA liquid crystals in vitro 
is vast and has been reviewed elsewhere (Livolant & Leforestier 1996). Durand et al (1992) 
performed careful X-ray diffraction experiments on slowly evaporating suspensions of 50 nm 
long DNA rods. The sequence of phase transitions was argued to be: isotropic - cholesteric 
- 2d hexagonal - 3d hexagonal - orthorhombic. The transitions are first order except for 
the one between the two hexagonal phases. In the 2d hexagonal phase, the rods are free 
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to move longitudinally but this freedom is gradually frozen out as the DNA concentration 
increases, so the 2d-3d transition is continuous. The solution only becomes orthorhombic at 
a very high concentration of 670 g/1 (the interaxial spacing H = 2.37 nm in the coexisting 
3d hexagonal phase). The helical repeat actually decreases quite dramatically from 3.5 nm 
to 3.0 nm within the sequence of positionally ordered phases, whereas the B conformation 
remains surprisingly intact with a constant rise = 2A = 0.336 nm. The buffer contained 0.25 
M electrolyte (either ammonium acetate or sodium chloride) so these data are pertinent to 
the packing of DNA under physiological conditions. 

The ionic-strength dependence of the cholesteric-hexagonal transition for DNA of length 
almost 50 nm was monitored by Kassapidou et al (1998), the hexagonal spacings at coexis- 
tence ranging from 3.7 nm at 0.10 M NaCl to 3.0 nm at 1.65 M. The orientational fluctuations 
measured by neutron scattering were close to 10 degrees. The first-order transition is visi- 
ble, in fact, in a test tube with a clear meniscus discernible between the two phases and the 
hexagonal phase exhibits the characteristic fanlike texture under a polarizing microscope 
(J.R.C. van der Maarel, private communication, 1998). 

On the other hand, Strey et al (2000) do not seem to witness the X-ray diffraction 
characteristic of hexagonal order, for concentrated Calf-Thymus DNA in 0.5 M NaCl, that 
has been wet-spun. One potential drawback is that the long DNA is quite polydisperse 
which could prevent a phase transition applicable to very long monodisperse DNA. Anyway, 
the azimuthal scattering has sixfold symmetry so that the phase is probably hexatic. If so, 
the decay in the orientational order should be algebraic (Nelson 2002) with an exponent 
depending on the conditions (like the ionic strength). This ought to be gleaned from the 
scattering curves (Lyuksyutov et al 1992). 

Spool models have been adopted as models of DNA packaged in bacteriophages (Harri- 
son 1973). X-ray diffraction patterns from phages have often been interpreted in terms of 
hexagonal packing of the DNA (North & Rich 1961; Earnshaw & Harrison 1977; Stroud et 
al 1981; Ronto et al 1988; Booy et al 1991; Cerritelli et al 1997; Bhella et al 2000). Cerritelli 
et al (1997) have adduced quite convincing evidence for the validity of the spool form for 
well- aligned bacteriophage T7. Earnshaw and Harrison (1977) pointed out that the inner 
region of the DNA spool should be more disordered. A theory of the nonuniform density 
of the DNA inside phages was set up recently (Odijk & Slok 2003) in which a cholesteric- 
hexagonal transition could be exhibited even under conditions of rather tight packing (note 
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the spacings quoted above). Physicists have analyzed the X-ray scattering by nematic liquid 
crystals of poly-7-benzyl glutamate in terms of density fluctuations (Ao et al 1991). These 
are coupled to director fluctuations by so-called continuity laws (de Gennes 1977, Selinger & 
Bruinsma 1991, Kamien et al 1992). The intensity spectrum as a function of the scattering 
vector has a characteristic bowtie shape. There is a typical length scale associated with these 
fluctuations which is determined by a trade-off of bending versus osmotic compressibility 
(de Gennes 1977). This turns out to be similar in size to the radius of the inner hole of 
the spool inside phages (bulk hexagonal DNA here being compared to phage DNA). Thus, 
it would be of some interest to reinvestigate the scattering by the inner region of the DNA 
spool from this point of view (Odijk & Slok 2003). 

The osmotic pressure of DNA solutions has also been readdressed recently (Strey et 
al 1999, Raspaud et al 2000). At very high concentrations where the phase is hexagonal 
or hexatic (but 1 shall assume it to be hexagonal for definiteness) , there are significant 
deviations from the Oosawa theory (eqs (17) and (18)). The approximate theory in terms 
of a hollow cylinder and thin boundary layers goes some way towards explaining the high 
values of n os (Strey et al 1999). For instance, at a DNA concentration of 1.2 monoM, eq (27) 
predicts 2.4 x 10 6 N/m 2 in good agreement with experiment (eq (31) does not apply since 
the separation is h — 0(a) and not h <^ a). At a concentration of 1.8 monoM, eqs (27) and 
(31) yield osmotic pressures 5.2 x 10 6 N/m 2 and 13 x 10 6 N/m 2 respectively, compared with 
the experimental 13 x 10 6 N/m 2 (Strey et al 1999). 

Till now, we have neglected configurational entropy. The free energy of a wormlike chain 
trapped in a tube of diameter less than the persistence length P was formulated long ago 
(Odijk 1983) and has been discussed many times since. For DNA confined within a hexagonal 
lattice, we may write on the basis of scaling arguments (Odijk 1983, Selinger & Bruinsma 
1991) 

fen(S) = j^tj-, (40) 

n ° s = 33/2^5/3^1/3 ( 41 ) 
Unfortunately, the coefficient c\ is not at all known accurately for a hexagonal phase. Setting 
the coefficient c± — 1 and using eqs (7), (12) and (14) lead to a semiquantitative explanation 
of the fairly sharp rise in the loading force as an elastic worm becomes densely packed 
in a cavity (see the simulations of Kindt et al (2001) and of Marenduzzo and Micheletti 
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(2003). The latter authors also find that there is a section of chain aligned within the 
cylindrical hole of the spool. Arsuaga et al (2002) put forward evidence for an inner region 
misaligned with the well-ordered outer spool under conditions of intermediate packaging 
but the effects of realistic electrostatics on these computations need to be awaited.) For 
comparison, the osmotic pressure of the entropic interaction is 5.3 x 10 5 N/m 2 at 1.8 monoM, 
an order of magnitude smaller than the pressures quoted above (h = 0.51 nm, P = 50 nm). 
Nevertheless, the nonlinear coupling between electrostatics and undulatory entropy probably 
needs looking into (see the work of De Vries (1997) on undulating charged surfaces in the 
salt-free Poisson-Boltzmann approximation). 

Loading forces and inner radii 

At high degrees of packing, the terms containing U in eqs (12), (13) and (14) are negligible 
so the relevant thermodynamic parameters reduce to V ~ n os ~ PksT /2SE 2 and \T\ = 
f+Sn os . The hole inside a spoollike region may have caps belonging to the wall of the capsid, 
a fact disregarded in section 2, but this effect gives only corrections of O [e 2 jV^ 3 j . We thus 
regain eq (10) of Odijk (1998) valid for a specific model, which itself is in good agreement 
with the DNA spacings measured inside bacteriophage T7 by Cerritelli at al (1997), Odijk 
(1998) employed the simple Oosawa model (eqs (17) and (18)). In recent papers, Tzlil et al 
(2002) and Purohit et al (2003) have used a different form: n os = ir [exp(H — H)/\ H — 1] 
with 7T = 0.12 fc B T/nm 3 , decay length \ H = 0.14 nm and H = 2.8 nm stemming from 
the osmotic stress measurements of Rau & Parsegian (1992) on DNA greatly perturbed by 
trivalent cobalt hexammine ions. Why this interaction, attractive in part, should generally 
prevail inside phages is not clear. Anyway, the spacings predicted by eqs (7) and (12) upon 
inserting this form, are virtually identical to those in table 1 of Tzlil et at (2002). The inner 
radius Ei is now 2.81 nm when the entire genome is packaged. 

The supposition of uniform density has been released in a recent analysis based on eq. (3) 
(Odijk & Slok 2003). In a density functional approach, the gradient terms are small because 
they exert themselves on a scale A or H only. The curvature stress is balanced against the 
chemical potential of inhomogeneous bulk DNA. The hole size is given by an expression 
similar to eq (12) but the tightly wound region surrounding the hole has an inhomogeneous 
density varying on a scale E { (Odijk & Slok 2003). Such a profile seems to be in qualitative 
agreement with experiments of Olson et al (2001) on the T4 phage. 
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The forces of encapsidation have been monitored in single-molecule experiments by Smith 
et al (2001). I shall try to interpret this force semiempirically. It seems a good Ansatz to 
set the osmotic pressure 7t os //cbT = p i = c 2 p = c 2 /AS within a certain interval of the DNA 
concentration in the absence of salt (see eqs (27) and (31) here, and fig 2 of Hansen et al 
(2001)). We then need the quantity f(S) at a certain ionic strength of the buffer. There 
is a Donnan equilibrium given by eq (17) and the free energy has been computed by Odijk 
& Slok (2003). We therefore obtain for the force \F\ = f(S) + Sn os = (c 2 k B T/A)[l - Inw] 
with w = 2c s /f>i = 2ASc s /c 2 <C 1. Smith et al (2001) measured a maximum loading force of 
50 pN at a buffer concentration c s of about 0.1 M. The inner volume of the phi29 capsid can 
be approximated by an ellipsoid of major axis equal to 25.5 nm and minor axis equal to 19.5 
nm (Tao et al 1998; P.J. Jardine, private communication, 2000). Hence we have Vq = 40600 
nm 3 and so the unit cell has dimensions S ~ 6.2 nm 2 and H ~ 2.7 nm for a genome length 
L max if we disregard the very small inner hole of size E i = 2.4 nm, computed from eq (12) 
upon a first iteration. We finally calculate c 2 = 0.73 and the osmotic pressure n os is about 
29 atm. This underestimates the theoretical and experimental bulk pressures a bit but the 
effect of magnesium ions in the force set-up has not been accounted for. 

However, there are difficulties in the interpretation of the force-length curve below the 
limit of full packing. If the inner hole remains fairly small, we may express the derivative as 

dL - L OS W 

using eqs (7) and (14). Accordingly, d\J-\/dL = ciksT j AL = 18 pN per L max if tt os = 
Cipo- This is not unreasonable if L ~ L max (see Smith et al 2001) but fails utterly when 
0.4 < L/L max < 0.9, for d\J-\/dL increases upon compaction in that case. It appears that 
we must have the osmotic pressure increasing as pi at least. What could be the cause of 
such a high effective power, granting that the measured force does not reflect dissipative 
losses? First, the crossover concentration demarcating the validity of the Oosawa model 
vis-a-vis a model with complete hexagonal detail is about 0.8 monoM judging from fig 2 of 
Hansen et al (2001) (L = 0.5L max at this concentration). Thus, counterions are gradually 
decondensing within the anomalous range of contour lengths. In addition, the Donnan 
effect is also slowly diminishing and the effect of the thin boundary layers is also becoming 
appreciable. An increase of the pressure stronger than linear is not surprising though this 
issue merits further quantitative study. (Parenthetically, a problem possibly related to this 
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could be the curvature in the rate-length curves determined by P.J. Jardine, opposite to 
the one measured by Smith et al (2001) (P.J. Jardine, private communications, 2000, 2001). 
He packaged DNA in phi29 phages in vitro by adding ATP to the solution. The reaction 
was effectively blocked with the help of gamma-S-ATP and the length of DNA injected was 
quantified as a function of time using a restriction enzyme.) 

The inner radii computed above at a high degree of packing are very small. An important 
problem that now arises is whether such tiny radii are viable, both physically and biolog- 
ically. We are here concerned with the tight curvature of long sections of solenoidal DNA 
independent of the base pair sequence and in the absence of bound proteins. Crick and Klug 
(1975) already quoted values of 3 to 5 nm for radii of curvature that could pose potential 
problems with regard to smooth bending. Sussman & Trifonov (1978) estimated a critical 
radius of kinking of about 2.5 nm by equating the elastic bending energy to the maximum 
enthalpy of stacking the base pairs assuming one kink per helical repeat. A more detailed 
treatment based on Van der Waals contacts led to a minimum acceptable radius of about 
5 nm. More recently, a 51 bp microcircle with a radius of curvature of 2.8 nm was devised 
by exploiting current computer modelling techniques (Tung & Soumpasis 1995). 42 bp mi- 
crocircles have actually been constructed in the laboratory (Wolters & Wittig 1989, 1990) 
but the base pair sequence is from a promotor. On the whole, a viable minimum radius 
would thus seem to be 3 nm or only very slightly smaller (A. Travers, private communica- 
tion, 2003). The inner radii Ei within phages that we have computed above are dangerously 
small. It is possible that the energy of curvature becomes anharmonic at small radii which 
would increase estimates of Ei somewhat. Any kinking would, of course, render invalid the 
models outlined here. 

DNA self-friction 

Jary and Sikorav (1999) investigated the monomolecular collapse of lambda DNA caused 
by the trivalent polyamine spermidine and the attendant rate of cyclization of the DNA. 
Fortunately, the collapse is very fast so it is possible to monitor the cyclization in the toroidal 
globular state quite accurately. The measured rates were between about 1 to 10 s _1 . This 
compares favourably with the reptation time from eq (37) (the viscosity of water t]q = 1 cP 
at room temperature, the contour length L — 17 jum and R t — 60 nm (Widom & Baldwin 
1983)). We conclude that the DNA apparently experiences friction from the "water" alone, 
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despite the fact that the typical spacing in such a globule is small (H ~ 2.8 nm). 

In a tightly packed phage, the dynamics of the DNA could be markedly different from that 
in globules, because of the much higher bending stress. In evaluating the frictional forces 
between two surfaces separated by a very thin layer of intervening fluid, one introduces a 
dimensionless group, i] u/tt os in this case, where u is the velocity tangential to the surface (see 
fig 5; Persson 2000). In DNA encapsidation within a phage, this velocity is so small (Smith et 
al 2001), that this quantity is minute so we ought te have boundary lubrication rather than 
hydrodynamic lubrication. The friction coefficient uo would be about 10 _1 (Persson 2000). 
In principle, there should then be a tendency for the water to be squeezed out between two 
adjacent DNA strands. Because there are only a few water molecules between the strands, 
layering transitions would be possible (Zilberman et al 2001) depending on the degree of 
commensurability of the two DNA backbones (Miiser et al 2001). However, eq (39) would 
now predict a minimum force more than 100 times larger than the 50 pN force observed in 
the phi29 phage. A possible resolution of this paradox is the very recent assessment of sliding 
friction between surfaces with intervening aqueous nanolayers (Raviv et al 2001, Raviv & 
Klein 2002). Surprisingly, the water still acts almost like a hydrodynamic bulk fluid, possibly 
because the hydrogen bonds are suppressed under conditions of sliding nanofriction. The 
friction coefficient could now be as low as 10~ 3 so the 50 pN force (Smith et al 2001) may 
be just enough to overcome the minimum force. It is also remarked that the free energy 
of compacting DNA into the phi29 phage is about 1 ksT per bp (Smith et al 2001) to be 
compared with the enthalpy of expulsion, 0.80 /c#T per bp, measured by calorimetry for 
bacteriophage T7 (Raman et al 2001; Serwer 2003), though the packing conditions in these 
two phages are not identical. Obviously, the problem of dissipative losses, if any, merits 
further attention. 

6. Concluding Remarks 

Clearly, our understanding of DNA packing in bacteriophages is incomplete. Although 
we may rationalize the high osmotic pressures and loading forces as being caused by the 
small separation between DNA helices, the tight inner radii of the spools are too close to 
kinking for comfort. At intermediate degrees of packing, we lack a precise theory of the 
electrostatic interaction in terms of the fraction of counterions condensed on the DNA. The 
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loading forces are very sensitive to detail in the interaction in this regime. An open problem 
is the friction at a high degree of compaction and its nature - is it hydrodynamic lubrication? 
Lastly, magnesium ions are often present in the buffers used. Their impact on the Donnan 
equlibrium needs to be investigated. 
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FIG. 1: A capsid packed with DNA down to an inner hole of volume Vq. 
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FIG. 2: The variational decrease in the size of the hole within the DNA spool. 
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FIG. 3: Three DNA cylinders within the hexagonal lattice. 
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FIG. 4: Cross section of the toroidal DNA globule. 
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FIG. 5: A standard arrangement to illustrate Coulomb friction. 
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